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Abstract 

We study the constraints that the operator product expansion imposes on large N c inspired 
QCD models for current- current correlators. We focus on the constraints obtained by going 
beyond the leading-order parton computation. We explicitly show that, assumed a given 
mass spectrum: linear Regge behavior in n (the principal quantum number) plus corrections 
in 1/n, we can obtain the logarithmic (and constant) behavior in n of the decay constants 
within a systematic expansion in 1/n. Our example shows that it is possible to have different 
large n behavior for the vector and pseudo-vector mass spectrum and yet comply with all 
the constraints from the operator product expansion. 
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1 Introduction 



The operator product expansion (OPE) has been used since long in order to gain information 
on the non-perturbative dynamics of the hadronic spectrum and decays [H [21 [3j HI El [6j 
[71 [HI [9]. In this article we revisit this problem. We want to obtain the constraints that 
the knowledge of the perturbative expansion in a s (Q 2 ) of the current-current correlators 
in the Euclidean regime poses on the relation between the decay constants and the mass 
spectrum for excitations with a large quantum number n (where n is the the quantum 
number of the bound state). We put special emphasis in going beyond the leading-order 
parton computation. We will work with a specific model for the hadronic spectrum. This is 
compulsory, since different spectral functions^ may yield the same OPE expression, yet we 
believe some aspects of our discussion may hold beyond the assumptions of our model. 

In order to have a well defined bound state it is crucial to consider the large N c approx- 
imation [10J. This ensures infinitely narrow resonances at arbitrarily large energies. We will 
consider to be in the large N c limit in what follows, as well as in the exact chiral (massless) 
limit. We will then set a specific model for the hadronic spectrum, valid for large values of 
n (we only need the behavior of the spectrum and decays for large n, we do not aim to get 
any information from perturbation theory for low values of n). This model will be based 
on the Regge behavior plus corrections in 1/n that will be included in a systematic way. 
The model is based on the assumption that the Regge behavior is a good description of the 
spectrum for large n (this can be explicitely seen in the 't Hooft model [UJ and it is also 
consistent with phenomenology). Given the 1/n corrections to the mass spectrum, the ex- 
pression of the correlator can also be written as a systematic expansion in 1/n, where higher 
powers in 1/n are equivalent to higher orders in 1/Q 2 in its OPE. By matching the OPE 
and hadronic expressions order by order in 1/Q 2 , we will be able to predict the logarithmic 
dependence on n of the decay constants (actually also the constant terms). This result can 
also be systematically organized within an expansion in 1/n together with an expansion in 
1/lnn. We will give explicit expressions up to order 1/n 2 and l/ln 3 n. We will also make 
some numerical estimates of the impact of these corrections. Finally, we would like to stress 
that we are able to introduce power corrections in 1/n to the Regge behavior and yet comply 
with the OPE. This is in contrast with Ref. [5], where, besides the Regge behavior, only 
exponentially suppressed terms are introduced (parametrically smaller than any finite power 
of 1/n for large n). This parameterization is however fine if considered as a fit not emanated 
from the large n limit. 



2 Correlators 

For definiteness, we will consider the vector-vector correlator but most of the discussion 
applies to any other current-current correlator (axial- vector, scalar, ....). 

nf (q) ee - g^q 2 )Tl v ( q ) ee i J d A xe^ x (vac\T {j£(x)J3)} \vac) , (1) 

1 In particular the one derived directly from perturbation theory, which we do not consider, since we will 
work in the large N c limit with infinitely narrow resonances. 
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where J v = (5/0/7'V/- In order to avoid divergences, we will consider the Adler function 

^).- Q >_n(Q*) = Q*/ o itjTTW -i M , (2) 

where Q 2 = —q 2 is the Euclidean momentum. 

Since we are working in the large N c limit, the spectrum consists of infinitely narrow 
resonances, and the Adler function can be written in the following way 

On the other hand, for large positive Q 2 , one may try to approximate the Adler function 
by its OPE, which reads 



(4) 



+ C{as ^ 2)) ' (3(a s (u))(vac\G 2 (u)\vac) + O (± 



where a^(Q 2 ) admits an analytic expansion in terms of a s (Q 2 ) (computed in the MS scheme), 

PW = -fo—4^-h[-4^-) +•••' ( 5 ) 
with p = 11/3JV C> pi = M/3N 2 , fa = 2857/54iV c 3 , and p2] 



2_ / _ 35 



C(a B (Q J )) = =- 1 - —Nc— 1 —-^- + ■■■). (6) 



3 Matching 

High excitations of the QCD spectrum are believed to satisfy linear Regge trajectories: 

M Vn 

lim — = constant. 

n— *oo n 

For generic current-current correlators, such behavior is consistent with perturbation 
theory in the Euclidean region at leading order in a s if the decay constants are taken to be 
"constants" , ie. independent of the principal quantum number n. 

The inclusion of subleading effects in a s can be incorporated into this model by changing 
the n dependence of the decay constants without changing the ansatz for the spectrum. The 
inclusion of these effects has consequences on subleading sum-rules and the relation with the 
non-perturbative condensates. 

Here we would like to go beyond the analysis at leading order in a s , as well as to consider 
power-like corrections in 1/n. We will consider that the large n expression for the mass 
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spectrum can be organized within a 1/n expansion in a systematic way starting from the 
asymptotic linear Regge behavior. In order to fix (and simplify) the problem we will assume 
that no Inn term appears in the mass spectrum^]. Therefore, we write the mass spectrum in 
the following way (for large n) 



M 2 (n) = B ( y s) n^ = B$>n + 4° } + + • • • , (7) 

s=-l 

where By - ^ are constants. We will usually denote M VLO (n) = By^n, M VNLO (n) = By^n + 

R (0) 
Cy. 



B { y ] and so on. To shorten the notation, we will denote By = By, By" 1 = Ay and By 



For the decay constants, we will have a double expansion in 1/n and 1/ Inn. 



F v( n ) = 2^ F vA n )— = F v.o( n ) + + — 

s=0 

where the coefficients F Vs (n) have a logarithmic dependence on n: 

oo ^ 



+ 



4>)=E4>)^- (9) 

r=0 

F 2 (n) 

As we did with the masses, we will define F VLO (n) = Fy (n), F VNLO (n) = F vo (n) + — , 
and so on. Note that in this case we also have an expansion in 1/ In n. 

We are now in position to start the computation. Our aim is to compare the hadronic 
and OPE expressions of the Adler function within an expansion in 1/Q 2 , but keeping the 
logarithms of Q. In order to do so we have to arrange the hadronic expression appropiately. 
Our strategy is to split the sum over hadronic resonances into two pieces, for n above or 
below some arbitrary but formally large n*, such that Aqcd^* *C Q. The sum up to n* 
can be analytically expanded in 1/Q 2 and will not generate InQ 2 terms (neither a constant 
term at leading order in 1/Q 2 ). For the sum from n* up to oo, we can use Eqs. (J7J) and (jHJ) 
and the Euler-MacLaurin formula to transform the sum in an integral plus corrections in 
1/Q 2 - Whereas the latter do not produce logarithms, the integral does. These logarithms of 
Q are generated by the large n behavior of the bound states and the introduction of powers 
of 1/n is equivalent (once introduced in the integral representation, and for large n) to the 
introduction of (logarithmically modulated) 1/Q 2 corrections in the OPE expression. 

Therefore, by using the Euler-MacLaurin formula, we write the Adler function in the 
following way (B 2 = 1/6, B 4 = —1/30, ...) 



"f 1 F y{n) r' F v (n) 

^ (Q 2 + M 2 (n)Y J (Qi + M 2 (n)y 



Q 2 F v( n *) | o2 fv , )k \B2k\ F 2 (n) 
2 (Q2 + M 2(n*)) 2-W j 



, (10) 

n=n* 



(2k)\dn^ 1 ) (Q 2 + M 2 (n)) 2 

2 This is a simplification. If one considers, for instance, the 't Hooft model [11], Inn terms do indeed 
appear. 
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where n* stands for the subtraction point we mentioned above, such that for n larger than 
n* one can use the asymptotic expressions d7j) and (|S]). This allows us to eliminate terms 
that vanish when n — > oo. Note that the last sum in Eq. fllOp is an asymptotic series, and 
in this sense the equality should be understood. 

Note also that for n below n*, we will not distinguish between LO, NLO, etc... in masses 
or decay constants, since for those states we will not assume that one can do an expansion 
in 1/n and use Eqs. (j7j) and (jSJ). 

Finally, note that the expressions we have for the masses and decay constants become 
more and more infrared singular as we go to higher and higher orders in the 1/n expansion. 
This is not a problem, since we always cut off the integral for n smaller than n* . Either way, 
the major problems would come from the decay constants, since, in the case of the mass, Q 2 
effectively acts as an infrared regulator. 



3.1 LO Matching 

We want to match the hadronic, Eq. fflUj) . and OPE, Eq. (j4]), expressions for the Adler 
function at the lowest order in 1/Q 2 . This means that we have to consider the lowest order 
expressions in 1/n for the masses and decay constants, i.e. F VLO {n) and My LO (n), since 
the corrections in 1/n give contributions suppressed by powers of 1/Q 2 - 

Only the first term in Eq. (fit)]) can generate logarithms or terms that are not suppressed 
by powers of 1/Q 2 - Therefore we obtain the following equality, 

This equation can be fulfilled by demanding that 

Fv,w(") 1 



I^W/*!-,^^'"')' (12) 

By using the perturbative expression for Imll^f r " (see [13]), we obtain 

F -- M = Bv It^ + i N - a - (nBv) + 243 i 28 7 ^ c(3) N ' a ' {nBv) 

f k 

346201 - 2904tt 2 - 324528 C(3) + 63360 C(5) Ar3 3 , _ , ^ . 4 . _ , x l 
+ 276487r3 — K^{nBv) + O {at(nB v )) j 

where a s (nBv) should actually be understood as a function of a s {By) and Inn. Therefore 



(13) 
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it is obvious that the above expression is resumming powers of a s (By) Inn: 



4 N r . 



3 16tt 2 



J2q} i+ 



2 1 + iA Nc ^xl Hn 



N 



aJB 



v 



47T 



(14) 



(2673 - 1936 C(3) - 408 In (l + §-N c a s {B v ) ln(n))) N 2 c a 2 s (B x 



N^(B X 



11 at a s (B v ) 
3 iVc 4tt 



ln(n))" 



(4tt) 3 52272vr(l + ^Ar c ^^4^1n(n)) 3 
+ 121tt (346201 - 2904tt 2 - 324528 C(3) + 63360 C(5)) 

-3672tt(2877 - 1936 C(3)) In ( 1 + ^-a s (B v ) m(n 



(An) 2 

350427 N c a s {B v )ln{n) 



12tt 
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+749088tt InM 1 + — a s (B v ) ln(n 

127T 



+ ^ (« s 4 (^)) 



Doing so we see that we are able to obtain the dependence of the decay constant in Inn 
(somewhat we are assuming that a s (By) is an small parameter, By ~ 1 GeV). 

We can also rewrite the decay constant as an expansion in 1/ Inn by using the equality 



Inn 



4?r + fa ln f a s (nB v 
(3 \a s {nB v ) ft \° 47r 





a s {nB\ 
4ir 



(15) 



where n = nBy / 'Aj^g. We then obtain 



Fyioin) 



B, 



4 Nr. 



3167T 2 

1 



ln 3 n 



/ 

46818 
161051 



9 1 
22 LT^ 

459 



In 2 n 



459 , , „ 144 /243 11 , , 

In Inn + C(3) 

1331 121 V 128 8 v ; 



42272605 

ln z lnn+ — (-2877 + 1936C(3))lnlnn+ 

322102 v sv ;; 2576816 



3vr 2 20283 C(3) + 360 C(5) 



22 



1331 



121 



O 



1 



In n 



(16) 



Note that the lowest contribution in 1/lnn for the decay constant, By^^s $2/ Q}> 
which, usually, is the only one considered, reproduces the leading-order partonic prediction 
for the Adler function. 

Note also that there is no problem with the Landau pole, even if the result is written in 
the form of Eq. fflBI) . since it holds only for n larger than an n* such that A^jg -C n*By (the 
integral has an infrared cutoff at n*). 

Finally, we remind that, strictly speaking, we can only fix the ratio between the decay 
constant and the derivative of the mass. We have fixed this ambiguity by arbitrarily imposing 
the n dependence of the mass spectrum. 



3.2 NLO matching 

We now want to obtain extra information on the decay constant by demanding the validity 
of the OPE at 0(1/ Q 2 ), in particular the absence of condensates at this order. We then 



6 



have to use the NLO expressions for M v {n) and F v {n). With the ansatz we are using for 
the mass at NLO, it is compulsory to introduce (logarithmically modulated) 1/n corrections 
to the decay constant if we want this constraint to hold. Note that it is possible to shift all 
the perturbative corrections to the decay constant. 

Imposing that the 1/Q 2 condensate vanishes produces the following sum rule: 



A A 1 

A— — A pt - — A pt + — 



n*-l 
n=0 



]k \B 2k \ d^-V 
Q 2 ^ (2k)\dn( 2k - 1 ) ±v 



F'(n) 



k=l 

oo 



Q 2 



d nF v,Lo( n ) 



dn 



+ 



F v (n*) 
2Q 2 

F v,i( n )/ n 
(Q 2 + M VjLO (n))< 



(17) 



+Q 5 



dn 



F VA (n)/n 



o (Qi + M V r (n)y 



0. 



This equality should hold independently of the value of n*, which formally should be taken 
large enough so that a s {n*By) <C 1, i.e. the limit A^jg- <C n*B v <C Q 2 ■ Again, the meaning 
of the asymptotic series appearing in Eq. ( flTj) should be taken with care. If we forget about 
this potential problem, only a few terms in Eq. ( TT7T) can generate In Q 2 terms, which should 
cancel at any order. Those are the first two and the last two terms. Actually, the next to last 
term does not generate logarithms, but it allows to regulate possible infrared divergences 
appearing in the calculation. Therefore, asking for the cancellation of the 1/Q 2 suppressed 
logarithmic terms produced by the first two and the last term in Eq. (flTl) fixes F vl . The 
non-logarithmic terms should also be cancelled but they cannot be fixed from perturbation 
theory. 

One can actually find an explicit solution to the above constraint for F V1 by performing 
some integration by parts. We obtain 



FyAn) 



n 



Ay d 

BV,d^ Iv ' 



F 2 Jn) 



A 



4 N c 



v 



3 16vr 2 



9 1 
22 hFft 



459 , w 2187 36C(3) 

(1 -2 In (Inn)) + 



1331 



484 



11 



ln 3 n 



2576816 



[-45794053 + 351384vr 2 + 41637552 C(3) - 7666560 C(5) 



-3672 In (Inn) (-3013 + 1936C(3) + 204 In (Inn))] + 



In n 



ln 5 n 



or in terms of a s {nBy) or a s (B\ 



FyAn) 



Ay^ 



4 N r 



n 



3167T 2 



11 , T 2 2r r, \ 2877- 1936 C(3) Ar3 3 . „ 
N 2 a 2 s (nB v ) 76 87r 3 N ^ nB ^ 



32tt 



11(376357 - 2904tt 2 - 344112 C(3) + 63360 C(5)) 
1105927T 4 



(19) 



N^at{nBy) + O (a 5 s (nB v )) 
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F vM . A 4 N c v 1 ( 11 1 2 a, 2 (i?y) 

^3167^^71 \ 2(1 + ^^)^)2 c (47r)2 I 

(2877 - 1936 C(3) - 408 In (l + ^ N c a s (B v ) ln(ra))) N3 a*(B v ) 



12(l + fN c ^±\n(n)) 3 ° (4tt) 3 

' N t a u B \) [-233618 N c a s (B v ) ln(n) 



4752vr(l + ^iV c ^?^ln(n)) 4 ° (4?r 
+ 121tt (376357 - 2904tt 2 - 344112 C(3) + 63360 C(5)) 

-3672tt(3013 - 1936 C(3)) In ( 1 + -^-a s {B v ) ln(n 



12tt 

+749088tt In 2 (\ + ^-a s (B v ) ln(n)^ + O (a 5 s (B v )) 

Note that besides the 1/n suppression, we also have an extra a 2 (nB v ) suppression. 

In principle one could think of the existence of 1/n x constant terms in the decay constant, 
i.e. without any associated logarithm. However, such terms produce ln(Q 2 ) contributions in 
the Euclidean regime that do not appear in the perturbative computation, so they can be 
ruled out. This appears to be a general statement since l/n m x constant for any m integer 
also produces logarithms. Note that in order to give meaning to these integrals it is implicit 
that the integral over n has an infrared cutoff at n*. Nevertheless, the logarithm does not 
appear to multiply powers of the infrared cutoff (as expected). 

Finally, we would like to mention that, besides the constraints coming from the logarith- 
mic related behavior of the OPE, there is also the constraint from its constant terms, which 
should sum up to zero. Nevertheless, for this constraint we cannot give a closed expression. 
This is due to the fact that the InQ 2 independent terms may receive contributions from 
any subleading order in the 1/n expansion of the masses and decay constants. The reason 
is that the decay constant at a given order in 1/n is obtained after performing some inte- 
gration by parts, which generate new (In Q 2 -independent) terms that can be Q 2 enhanced. 
This statement is general and also applies to any subleading power in the 1/Q 2 matching 
computation. 



3.3 NNLO matching 

We now consider expressions for the mass and decay constants at NNLO. For the first time 
we have to consider condensates. Simplifying terms that do not produce logs, we obtain the 
following equation, 



35 a s (Q 2 ) P(a s {v)){vac\G 2 (v)\vac) 



121 4tt 

12 



Q 4 



(21) 



= Q 2 



dn 



B . (Q 2 + B v nf 



n 2 



1 d 

By dn 



C v F 2 fi (n)_ + A v F 2 ^n) 



n 



2n 
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where = stands for the fact that we can only predict the InQ 2 dependence. Constant terms 
are not fixed by this relation. 

In order to get a more closed expression is convenient to use the following equality, 



1 



a s (Q 2 ) = Q 2 



dn 



(Q 2 + B v n) 2 B v n 2 8vr 



^° 2/ D \ 

a s [riBy) 



(22) 



valid up to terms that do not produce logarithms or those that are subleading. 
We get then 



V,2 



n 



-Cv~: 



4 N r 



3 16tt 2 



J2Q 2 f\ 1 + t Ncas{nBl 

f v 



287- 176 C(3) _ 

128vr 2 6Att 2 B v C v 
+0 (a 3 s (nB v ))} . 



35 (3 {a s {v))(vac\G 2 {u)\vac) 
88 B V C V N 2 



N 2 a 2 (nB v ) (23) 



Note that in this case we only consider up to 0(a 2 (nBy)) corrections, since higher order 
loops are unknown. The accuracy is set by our knowledge of the matching coefficient of the 
gluon condensate. Note also that F V2 (n) does not have a s suppression. Therefore, for low n, 
this contribution could be practically of the same size than, formally, more important terms. 



4 Axial versus vector correlators 

The above discussion has been performed for the vector-vector correlator Adler function. 
It goes without saying that we could perform a similar analysis with axial-vector currents, 
since the perturbative expansions for both correlators are equal. Here it comes an important 
observation. We could change the coefficients for the mass spectrum Ba ^ By, A a ^ Ay, 
• • • , yet we would obtain the same expression for the OPE (at the order we are working, the 
first chiral breaking related effects are 0(1/Q e )). Therefore, we conclude that the OPE does 
not fix B A = By as it is sometimes claimed in the literature [JJ [5ffi Our computation gives 
a specific counter example. Moreover, it is nice to see what the role played by Ba and By 
is in our case. When one goes to the Euclidean regime, Ba and By become renormalization 
factorization scales and, obviously, the physical result does not depend on them (for large 
Q 2 in the Euclidean). On the other hand, it is evident that having different constants: Ba, 
By, . . . produces different physical predictions for the masses and decay widhts for vector or 
axial-vector channels. The point to be emphasized is that Ba = By cannot come from an 
OPE analysis alone. This point has already been stressed in Refs. jU E], what we think is 
novel in our analysis is that we have seen that the inclusion of corrections in a s does not 
affect that conclusion, and that Ba and By play the role of the renormalization scale in the 
analogous perturbative analysis in the Euclidean regime, and are therefore unobservable. 
Finally, we cannot avoid mentioning the analysis of Ref. [2] where, using AdS/CFT, they 
explicitly find Regge behavior with different slopes for vector and axial-vector channels. 

3 Another issue, on which we do not enter, is whether some other kind of arguments (relying on the specific 
model used), like semiclassical arguments, may fix those parameters to be equal. 
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In any case, even though the constants that characterize the spectrum can be different 
for the vector and axial- vector channel, they have to yield the same expressions for the OPE 
when combined with the decay constants. This produces some relations that we list in what 
follows. We first define t = Byn = B^n' and take n and n' as continuous variables. We then 
obtain the following equalities 



By B A 
^2 



- 

TV 



^-lmU p v ert (t) = f (t), 



1 F v ^n) _ 1 n,i(n')_d 



A V B 



n 



A A B 



ADA 



TV 



(24) 
(25) 



where 



1 

By' 



FvM 



1 

B~ A 
1 



1 d f C v F Vfl {n) | A v F Vtl (n) 
By dn \ n 2n 



TV 



r2 



/ C A Fj ( n') AaFI^u') 
Ba dn 1 \ n' 



+ 



2n' 



= ^(3(a s (u))(vac\G 2 (u)\vac)f 1 (t), 



35 (3 a 2 s (t) 
T2lT(47r) 2 



+ 



(26) 
(27) 



5 Numerical Analysis 





n = 1 


n = 2 


n = 3 


n = 4 


M p (l) 781.3(775.5 ±0.4) 1440.2(1459 ± 11) 


1891.8(1870 ±20) 2257(2265 ± 40) 


M p (II) 


771.5(775.5 ±0.4) 


1471.7(1459 ± 11) 1855(1870 ±20) 


2154.8(2149 ± 17) 


M ai 


1235.6(1230 ±40) 


1621.7(1647 ±22) 1962(1930^) 


2257.8(2270+^) 


F v (l) [ 156(156 ±1) 


155 154 


153 


F V (1T) 185(156 ±1) 


147 


139 


135 


F A 123(122 ±24) 


137 


139 


139 



Table 1: We give the experimental values of the masses (in MeV) and electromagnetic decay 
constants (when available) for vector and axial vector particles (within parenthesis) , compared 
with the values obtained from the fit. For the vector states we consider two possible Regge 
trajectories that we label I and II respectively. We take a s (lGeV) = 0.5 and (3{G 2 ) = 
-(352 MeV) 4 . 

We restrict ourselves to the SU(2) case (non-strange sector) and study the vector and 
axial- vector channels. We would like here to assess the importance of including perturbative 
corrections to a standard analysis based on the OPE. We do not aim to perform a full fledged 
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analysis, but only to see the importance of the corrections. In table [TJ we give the values 
of the masses and decay constants. In Figured] we show the changes in both F VjLO (I) and 
Fy t Lo{N) as we include higher orders in the expansion in a s , and in Figure [2] the changes 
in the full i*V(I) an d Fv(U) as we include higher orders in 1/n. In figure [3] we show the 
same plots for the axial- vector case. We take the experimental values from Ref. [15J . In 
principle there are more states in the particle data book, in particular in the vector channel. 
Nevertheless, it is not clear whether they belong to the same Regge trajectory or whether 
they belong to some daughter one, see, for instance, the discussion in Ref. [5]. For the time 
being we will disregard the study of other possible (vector) Regge trajectories and restrict 
the analysis to a single trajectory. We will consider the two possibilities listed in Table [TJ 
Our choice of states for the set (I) is motivated by the discussion of Ref. [16] on the possible 
formation of multiplets in the case of chiral symmetry restoration. The set (II) is based on 
the assignment of states made in Ref. [5] (based on the existence of S and D-w&ve daughter 
trajectories) and in particular on the analysis of Ref. [17], where the state 2265 is argued to 
belong to the D-w&ve Regge trajectorjo 

In order to fix the parameters of the mass spectrum we use the experimental values of 
the masses we list in the table. We obtain the values: 

B v (l) = 1.525 x 10 6 MeV 2 , A v (l) = -1.038 x 10 6 MeV 2 , C v (l) = 0.123 x 10 6 MeV 2 , 
B v (Ll) = 1.128 x 10 6 MeV 2 , Ay (II) = 0.353 x 10 6 MeV 2 , C V (II) = -0.885 x 10 6 MeV 2 , 
B A = 1.278 x 10 6 MeV 2 , A A = -0.100 x 10 6 MeV 2 , C A = 0.349 x 10 6 MeV 2 . 

(28) 

We should mention that the values obtained for these parameters are not very stable 
under the change of number of data points, except for By and B A , which are roughly stable, 
although with quite sizeable uncertanties. For the subleading terms A and C, their values 
are basically random with the fit. We roughly find By ~ B A within the uncertainties. The 
n dependence of the axial and vector (model II) decay constants is small but sizeable (and it 
goes in the right direction for low n). The 1/n corrections are always corrections compared 
with the leading order terms. Nevertheless, the 1/n 2 correction is much larger than the 1/n 
one for the range of values of n that we explore. This appears to be due to the a 2 /(4-7r) 2 
suppression of the 1/n term, as well as to the difference in size between the constants A 

4 We also thank S. Afonin for discussions on this point. 
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Figure 2: In this plot we show Fy(V) and Fy (11) at different orders in the 1/n expansion. 




Figure 3: In this plot we show Fa^lo an d Fa at different orders in a s and in the 1/n 
expansion, respectively. 



and C. This is so for the axial and vector (model II) decay constants. Nevertheless, for the 
vector (model I) decay constants the n dependence appears to be quite small also at NNLO. 
This appears to be due to the small value of the coefficient Cy(l). The gluon condensate 
contribution is a small correction to the total NNLO term. Either way, our predictions 
compare favorably with experiment when this comparison is possible. 

We should keep in mind that these results have been obtained for a specific model, so 
we are testing the impact of the perturbative corrections for this specific model. On the 
other hand, if one believes that the large n behavior of the spectrum is dictated by the 
Regge behavior and that the corrections can be obtained as an expansion in 1/n, the set 
up is general. The only ambiguity comes from where the logarithms should be introduced 
(masses or decays). At this respect it is worth mentioning that, as a matter of principle, this 
ambiguity could be fixed if enough experimental information were available for the masses 
and decays. 

6 Conclusions 

We have studied the constraints that the OPE imposes on large N c inspired QCD models for 
current-current correlators. We have focused on the constraints obtained by going beyond the 
leading-order parton computation. We have explicitly showed that, assumed a given mass 
spectrum (Regge plus corrections in 1/n), we can obtain the logarithmic (and constant) 
behavior in n of the decay constants within a systematic expansion in 1/n. More than that, 
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power-like 1/n corrections can only be incorporated in the analysis if full consideration to the 
perturbative corrections in the Euclidean regime is made. This is due to the fact that these 
type of contributions produce logarithms of Q in the Euclidean (this is one of the reasons 
why this sort of corrections are not usually considered in quark-hadron duality analysis). On 
the other hand, the existence of Inn in the decay constants may point to the existence of 
two scales in the problem, Aqcd and uAqcd, in the Minkowski regime. 

We have also performed some numerical estimates of the importance of these corrections. 
The n dependence of the decay constants is small but sizeable for the axial and vector (model 
II) channel, for the vector (model I) one this dependence is small. On the other hand the 
uncertainties of the calculation are large. Either way, our predictions compare favorably 
with experiment when this comparison is possible. 

Our example shows that it is possible to have different large n behavior for the vector 
and pseudo-vector mass spectrum and yet comply with all the constraints from the OPE. 

An important caveat of our analysis is that we have not considered what the effect of 
renormalons could be. We have focused on the effect of low orders in perturbation theory 
to the decay constants. It would be interesting to see whether the knowledge of the higher 
order behavior of perturbation theory may give some extra constraints on the values of these 
constants and the mass spectrum. At this respect we have to say that we have obtained 
approximated expressions for the decay constants as an expansion in a s (nBy), with just 
the low order contributions in a s . It is quite likely that this expansion is asymptotic and 
that different orders in 1/n are related in a similar way to the one found in the renormalon 
analysis for the OPE expansion for different orders in 1/Q 2 . Therefore, the results obtained 
for the 1/n corrections could be affected as well by the asymptotic behavior of the 1/lnn 
expansion in the leading-order term. This is obviously related with renormalons. We expect 
to come back to this issue in the future. 
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